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ABSTRACT 

In this paper we report on Chandra observations of the starburst galaxy 
NGC 922. NGC 922 is a drop-through ring galaxy with an expanding ring of 
star formation, similar in many respects to the Cartwheel galaxy. The Cartwheel 
galaxy is famous for hosting 12 ULX, most of which are in the star forming 



-2- 



ring. This is the largest number of ULX seen in a single system, and has led to 
speculation that the low metallicity of the Cartwheel (0.3 Z ) may optimize the 
conditions for ULX formation. In contrast, NGC 922 has metallicity near solar. 
The Chandra observations reveal a population of bright X-ray sources, including 
7 ULX. The number of ULX in NGC 922 and the Cartwheel scales with the 
star formation rate: we do not find any evidence for an excess of sources in the 
Cartwheel. Simulations of the binary population in these galaxies suggest that 
the ULX population in both systems is dominated by systems with strong wind 
accretion from supergiant donors onto direct-collapse BHs. The simulations cor- 
rectly predict the ratio of the number of sources in NGC 922 and the Cartwheel. 
Thus it would appear that the the metallicity of the Cartwheel is not low enough 
to see a difference in the ULX population compared to NGC 922. 

Subject headings: galaxies: individual (NGC 922) - X-rays: galaxies 



Introduction 



Chandra images of starburst galaxies are spe ctacular, showin g a multitude of bright 
point sources and plumes of X-ray emitting gas ( jFabbianol 120061 ) . The brightest of the 
discrete sources (L x > 10 39 ergs s _1 ) are known as Ultra-Luminous X-ray Sources (ULX), 
and have attracted considerable attention in recent years because they have broad-band X- 
ray l uminosities many times the Eddington limit for a neutron star or stellar mass black 
hole (IZezas. Fabbiano. Rots, fc Murray! 120021 e.g.). There is a growing concensus that most 
ULX can be explained as high (> 10 M ) stell ar mass black h oles radiating in excess of the 
Eddington limit, possibly with mild beaming ([Roberts! 120071 see) and references therein). 
The most massive stellar black holes are thought to form in low metallicity environments, 
where the stellar winds of the progenitor are relatively weak and do not carry of f such a 
large fraction of the initial mass of the star ( IBelczynski. Sadowski. fc Rasiol 120041 ) . Thus 
ULX should form preferen tially in low meta llicity environments. Observational evidence for 
this is seen in some ULX ( ISoria et al.ll2005[ e.g.) 



NGG 922 is a peculiar galaxy with a distinctive C-shaped ring of Ha e mission (see 



Figur e [T]) . It has been described as a dust-obscured g rand design spiral galaxy ({Block et al. 



20011 ). However, a recent paper by lWong et al.l (120061 ) presents compelling evidence that the 
morphology could be explained as a result of an off-center collision between a spiral and its 
dwarf companion. This model successfully explains the C-shaped morphology. The intruder 
galaxy is identified to be a dwarf galaxy approximately 8 arcminutes away, apparently con- 
nected to NGC 922 by a plume of tidal debris. The starburst in NGC 922 is still on-going 
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as evide nced by the luminou s Ha emission and a population of very young, massive star 
clusters (jPellerin et al.ll2010[ ). NGC 922 is similar to the Cartwheel in that both of them 
are the remnants of drop through collisions, and both of them have high star formation 
rates. The C artwheel galaxy is famous for hostin g 12 ULX, most of which are in the st ar 
forming ring ( jWolter fc Trinchieril (12004 ). see also iGao. Wang. Appleton. fc Lucas! (120031 )). 
This is the largest number of ULX seen in a single system, and has led to speculation that 
the l ow metallicity of the Cartwheel (0 .3 Z Q ) may optimize the conditions for ULX forma- 



tion dMapel 



( IWong et al 



i. Colpi. fc Zampieril 120091 ). In contrast, NGC 922 has metallicity near solar 



2006). We are currently engaged in a Chandra and HST study of the X-ray 



binary population of NGC 922. In this paper, we compare the ULX population of NGC 922 
and the Cartwheel. Our primary goal is to determine whether there is an excess of ULX 
(normalized to the star formation rate) in the Cartwheel compared to NGC 922 that can be 
attributed to the difference in metallicity. In addition, we use population synthesis models 
to estimate the expected number of ULX in both systems. In a future paper we use HST 
broad-band and Ha images to study the relationship between the X-ray sources and star 
clusters in NGC 922. 



2. Data Analysis 

Two images of NGC 922 were obtained in 2009. The first was a 30ks exposure designed 
to detect X-ray point sources down to a luminosity of 7xl0 38 ergs s _1 and the second a 
lOks exposure, scheduled a few months later, to detect variability in the brightest sources. 
Observational details are given in Table [2j Data were reduced using the software package 
CIAO (ver. 3.4) with the appropriate calibration products. We created exposure maps using 
the task mere;e_ and we determined the positions of the point sources using the CIAO 
task wavdetectj^] A smoothed X-ray image of NGC 922 is shown in Figure [2j The X-ray 
point sources are marked. Counts for the point sources were extracted from regions shown 
in Figure EJ Background counts for each point source were extracted from an annular region 
around each source (these regions are not shown in Figure [2]). The counts for each source 
were corrected for the PSF fraction (Point Spread Function fraction, i.e. the fraction of 
counts outside the extraction region) by using a PSF with energy 1 keV. The background 
subtracted, PSF corrected count rate for each source at both epochs is shown in Table [2j 

We extracted X-ray spectra and standard response files (ARF and RMF) for each 



1 See http://cxc.harvard.edu/ciao /threads / merge.all 

2 See http:/ /cxc. harvard. edu/ciao3.0/download/doc/dctcct_html_manual/Manual. html. 
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source. Fluxes and luminosities for the deep March observation were constructed in the 
0.3-8 keV band, assuming a ga l actic absorption A^=1.6xl0 20 cm -2 . Using the results of 



Swartz. Ghosh. Tennant. &: Wul (120041 ). we adopt an intrinsic source spectrum with photon 



index T = 1.7 to calculate luminosities. Fluxes were then calculated in XSPEC by rescaling 
the normalization of the power law model to the observed cou nt rates. Luminosities were 
calculated assuming a distance of 48 Mpc (jOlivares et al.ll2010l ). 



3. The X-ray Source Population of NGC 922 



by L 



NGC 922 has a population of bright X-ray point sources, including 9 ULX, defined 
x > 10 39 ergs s _1 , in the 0.3-8.0 keV band (there are 7 ULX in the 2-10 keV band). 
Count rates for both the March and October observations are shown in Table [2j Many 
of the sources are variable. Figure [1] shows a Ha image of NGC 922 with positions of 
the X-ray sources superimposed. The X-ray sources are clearly associated with recent star 
formation, and several are coincident with the star forming ring. The close association of 
bright X-ray sources with recent star formation has been observed in many starburst galaxie s 
(IWolter fc Trinchierill2004j ; IZezas et al.ll2007l ; iRappaport. Levine. Pooley. fc Steinhornll2010l ) . 
It is generally accepted that the X-ray source pop ulation at high luminosities (L x > 10 37 



ergs s 1 ) is dominated by accreting X-ray binaries (lFabbianoll2006l ) 



3.1. The Brightest Source in NGC 922 



NGC 922 X-4 has an 0.3-8 keV luminosity L r = 2xl0 40 ergs s" 1 , placing it at the 



upper end of the X-ray binary lum inosity function ( IGrimm. Gilfanov. &: Sunyaevl 12003 



Swartz. Ghosh. Tennant. &: Wull2004l ). The count rate declined by a factor of 2 in the sec- 
ond observation. We fit the 0.3-10 keV spectrum of NGC 922 with two models - a simple 



ObsID 


Date 


Exposure Time 


Chip 


Data Mode 


Exposure Mode 






ks 








10563 


2009-03-05 


29.737 


S3 


TIMED EXPOSURE 


VFAINT 


10564 


2009-10-02 


10.020 


S3 


TIMED EXPOSURE 


VFAINT 



Table 1: Observing configuration for two Chandra observations of NGC 922 
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Net Count Rate F x (0.3-8 keV) L x (0.3-8 keV) 
Source RA Dec March 09 Oct 09 1(T 15 ergs cm" 2 s -1 10 39 ergs s" 1 

J2000 counts ks" 1 March 09 March 09 



1 


02 


25 


06.84 


-24 


46 


50.05 


0.92 


0.41 


6.77± 1.25 


1.86 


2 


02 


25 


06.50 


-24 


47 


17.19 


0.30 


<0.1 


2.49 ± 0.78 


0.68 


3 


02 


25 


05.76 


-24 


47 


01.78 


1.57 


2.04 


14.04 ± 2.0 


3.87 


4 


02 


25 


05.63 


-24 


47 


52.58 


11.47 


4.78 


87.80 ± 4.64 


24.20 


5 


02 


25 


05.28 


-24 


47 


57.41 


0.49 


1.06 


4.30± 1.05 


1.19 


6 


02 


25 


05.21 


-24 


47 


05.79 


0.44 


0.80 


3.74 ± 1.04 


1.02 


7 


02 


25 


05.05 


-24 


46 


53.98 


0.21 


<0.1 


1.89± 0.79 


0.52 


8 


02 


25 


04.92 


-24 


47 


59.36 


1.49 


1.00 


11.27 ± 1.65 


3.10 


9 


02 


25 


04.89 


-24 


47 


10.66 


0.44 


0.29 


3.57 ± 0.95 


0.98 


10 


02 


25 


04.02 


-24 


47 


02.58 


1.46 


1.16 


12.11 ± 1.77 


3.33 


11 


02 


25 


03.82 


-24 


47 


19.08 


1.68 


3.28 


14.20 ±1.89 


3.91 


12 


02 


25 


03.74 


-24 


47 


57.45 


4.08 


3.49 


31.49 ± 2.81 


8.67 


13 


02 


25 


03.33 


-24 


47 


45.70 


0.42 


0.49 


3.23 ± 0.92 


0.88 


14 


02 


25 


03.32 


-24 


47 


29.45 


0.26 


0.21 


2.35 ± 0.77 


0.64 



Table 2: NGC 922 detected X-ray sources 



-6 - 



20 



power law and a multicolored disk blackbody. In both galactic of 1.62x10 

cm -2 was assumed and fixed in the fit, and another absorption component intrinsic to the 
source was allowed to vary. The resulting fit parameters are shown in Table 13.11 The 
luminosity, spectrum and variability of NGC 922 X-4 are typical of ULX that appear to 



have hard X-ray spectra at the highest acc retion rates ([Gladstone. Roberts, fe Pond 12009 



Berg. 



lea. Weaver. Colbert, fe Robertsll2008h . Other examples include I C 342 X-l, NGC 1313 



X-l( 



Gladstone. Roberts, fe Dondl2009h. NGC 1365 flSoria et all 200% and Holmberg IX X-l 



(IKaaret fe Fengj|2009i ; IVierdayanti. Done. Roberts, fe Mineshigdl2010l ). This high flux/hard 



spectrum behavior is interesting because it appears to be different to what is generally ob- 
served in Galactic black hole binaries, where soft thermal emission from a disk dominates 
higher accret ion states. The "high/hard" state as been provisionally called the " ultralumi- 
nous" state ([Gladstone. Roberts, fe Dondl2009l ). It may correspond to a super-Eddington 
accretion state where an out-flowing wind fro m the inner accretion disk leads to a n increased 
optical depth and lower coronal temperature (IGladstone. Roberts, fe Dondl2009i ). 



3.2. Comparison with the Cartwheel 



Here we compare the ULX population seen in NGC 922 with that of the Cartwheel. The 
Cartwheel has a metallicity in t he range Z~ 0.18 — 0.3Z Q in the ring (IFosbury fe Hawarden 
1977l ; IVorobyov fe Bizvaevll2001). Th is is considerably lower than the metallicity of NGC 922, 



which is ~0.75 Z Q (jWong et al. 



20061 ). Figure [3] shows the cumulative luminosity function for 



sources in NGC 922 and the Cartwheel. Luminosities for sources in the Cartwheel were taken 



from IWolter fe Trinchieril (120041 ) . In order to check that there were no s ystematic differences 
i n flux es derived using our method (described in [2]) and that used by IWolter fe Trinchieri 
(120041 ) we repeat ed our analysis using the ar chival Cartwheel data. We derived fluxes essen- 
tially identical to lWolter fe Trinchieril ( 120041 ) . We no te that luminosites for X-r ay sources in 
the Cartwheel are calculated in the 2-10 keV band (IWolter fe Trinchieril 120041 ). The X-ray 
luminosities of sources in NGC 922 have been plotted in the 2-10 keV band in order to allow 
for direct comparison with the Cartwheel. This results in slightly lower luminosities than 
given in Table H There are 7 ULX in NGC 922 in the 2-10 keV band. 

We fit a si ngle power law to the unbinned , diffe rential, XLFs using a maximum likeli- 
hood statistic (jCrawford. Jauncey. fe Murdoch! 1 19701 ) . The goodness-of-fit (GOF) estimate 
is performed by simulating a luminosity distribution with the best-fit slope. One million it- 
erations were performed. If the data are well fitted by a single power law, the GOF statistic 
approaches 1.0. The contribu tion of cosmic backgr ound sources to the XLF was evaluated 
using the log N-log S curves of iGiacconi et al.l (120011 ). Details of the fitting method are given 
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in iKilgard et al.l (120051 ) . Table I3~2l gives the results. 



Table 13.21 shows that the XLF slopes of high luminosity X-ray sources in the Cartwheel 
and NGC 922 are identical within the uncertainties. The slopes are consistent w i th th e 
"Universal" High Mass X-ray Binary XLF described by lGrimm. Gilfanov. fc Sunyaevl ( 120031 ). 
Assuming th at the star formation rates of the Car twheel and NGC 922 are 18 and 8 M Q /yr 
respectively (IMayya et al.ll2005l ; IWong et al.ll2006l ). we find that the number of ULX sources 
relative to the star formation rate is 0.94±0.3 ULX per solar mass of star formation for 
NGC 922, and 0.67±0.2 ULX per solar mass of star formation for the Cartwheel. Here we 
use the number of ULX in the 2-10 keV band. These ratios are dependant on the assumed 
star formation rates. The star formation rates quoted for NGC 922 and the Cartwheel were 
derived using extinction-corrected Ha luminosities. There is no evidence for widespread 
obscured star formation in either galaxy, so these values are likely to be close to the correct 
values. We conclude there is no evidence for an excess of high luminosity sources (normalized 
to the star formation rate) in the lower-metallicity Cartwheel unless the SFR of either galaxy 
is wrong by a factor of 2 or more. 



4. Population Synthesis Models 

The similar ULX/SFR observed in the Cartwheel and NGC 922 is, in fac t, expected in 



popul ation synthesis models of ULX formation in low metallicity environments. iLinden et al. 



(120101 ) studied the mechanisms for high mass X-Ray binary formation in starburst environ- 
ments, finding the ULX population to be dominated by systems moving through pathways 
powered by strong wind accretion from supergiant donors onto massive black holes formed 
in failed supernovae at both Z = 0.75 Z and Z = 0.28 Z . 



1 20081 ) simulations using 
§ we calculate a popula- 



Utilizing simulations from the Star Track code (iBelczynski et al. 
a SFR of 18 M yr^ 1 for the Cartwheel and 8 M yr" 1 for NGC 922 [ 
tion of 103 ULX for NGC 922 and 238 ULX for the Cartwheel. We note two results: (1) our 
default simulations overpredict the absolute number of observed systems by approximately 
an order of magnitude, and (2) our models similarly predict no metallicity dependence in 
the ULX/SFR ratio in the regime between 0.75 Z and 0.3 Z investigated in this study. 

Predictions regarding the absolute number of bright X-ray Binaries are particularly dif- 
ficult for population synthesis models. We note several obvious parameter space choices 



3 The continuous star formation period is simulated to last for 100 Myr. However, a negligible number of 
ULX are produced after 20 Myr, rendering the exact end of the simulation unimportant. 
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Parameter Powerlaw Multicolored disk 

N H ,Gai (cm 2 ) 1.62xl0 20 1.62xl0 20 

(cm 2 ) 2AtH x 10 21 3.4+^ x 10 20 

7=1.9^| T m (keV) = 1.3±0.3 

normalization 2.3±0.95xi(T 5 1.36±1.01xl0 -3 



X 2 /DOF 1.15 (20.79/18) 1.15 (20.7/18) 



Table 3: Best fitting parameters for source X-4. Models are phabs*phabs*pow and 
phabs*phabs*diskbb. 



Galaxy 


7 




GOF 


Nfit 


N 39 






ergs s _1 








Cartwheel 


0.71±0.120 


lxlO 39 


0.81 


12.0 


11.1 


NGC 922 


0.76±0.22 


5xl0 38 


0.71 


14 


7.5 



Table 4: Power law fits to the XLF. 7 is the fitted slope to the XLF correcting for background 
sources. L m j„ is the minimum luminosity included in the fit. GOF is the goodness-of-fit 
parameter, Nfu the number of sources included in the fit and N39 the predicted number of 
ULX (i.e. sources with L x > 10 39 ergs s _1 ) 
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which would allow a match between our simulations and observations. Because the super- 
giant donors have extremely large radii, the ULX formed through this pathway must be 
weakly bound - implying that small natal kicks imparted to massive BHs (e.g. a Maxwellian 
distribut ion with dispersion velocity 26.5 km s -1 ) can disrupt as many as 90% of ULX pro- 
genitors (ILinden et al.ll2010l ). Additionally, our definition of ULX assumes an average X-Ray 



luminosity of 1 x 10 erg s , however eccentricities may make these systems dimmer for 
a substantial portion of their orbit. Another plausible explanation concerns the treatment 
of super-Eddington accretion. In our default simulations we allow black hole accretors to 
power ULX at up to ten times the Eddington limit, so long as the donor winds are able to 
supply enough material to the black hole accretor. Our models do not describe the specific 
mechanisms which allow this. Severa l mechanisms in volving dyna mic instabilities h ave been 
postulated - and include leaky-disk ( IBegelmanl |2002| ) or beaming (IKing et al.ll200lf ) models 
which may imply a duty cycle at the 10% level. 

The ratio of ULX formation rates between NGC 922 and the Cartwheel places more 
meaningful constraints on the models than absolute numbers, as there are fewer parame- 
ters which can fine tune t he metallicity depen dence. In this case, we find models using 
the default assumptions in ILinden et al.l (120101 ) to correctly predict insignificant metallic- 
ity dependence in the regime from Z = 0.75 Z Q to Z = 0.18 Z . We emphasize that this 
relation falls naturally out of our simulations, as the ULX population at both metallicities 
is created through the same evolutionary channels. We note that if the metallicity of the 
Cartwheel is significantly lower than assumed here, then this assumption would no longer 
hold, and the Cartwheel would begin to be dominated by a second population powered by 
Roche Lobe overflow of main sequence stars onto fallback BH accretors (this model has been 
suggested as a mechanism to explain t he properties of the brightest source in the Cartwheel 
( iPizzolato. Wolter. fc Trinchieril |2010| )). Thus the consistency of the ULX/SFR ratio pro- 
vides some preference for a higher modeled metallicity within the Cartwheel. Guided by this 
consistency in the ULX/SFR in the observed metallicity range, we note several other robust 
results. First, the supergiant donors which power ULX activity die off quickly after star 
formation, and we find ~95% of luminous ULXs to be younger than 20 Myr. Secondly, the 
direct collapse pathway imparts almost no systematic velocity to the resultant X-ray binaries 
, and thus ULX are constrained to lie very close to recent star formation regions. However, 
these systems may gain significant spatial velocities if they are forme d within dense stella r 



environments where three body interactions become non- negligible (jMapelli et al. 1 1201 ih . 



Finally, we caution that the significant overabundance of ULX sources predicted by Star- 
Track models implies an additional physical mechanism which must eliminate a large number 
of possible ULX systems. We are unable, in this study, to constrain this mechanism from 
introducing a metallicity dependence, which would create a mismatch between observations 
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and our simulated population. 



5. Summary and Conclusions 



NGC 922 is a drop-through ring galaxy ([Wong et al.ll2006l ) with an expanding ring of 
star formation, similar in many respects to the Cartwheel galaxy. New Chandra observa- 
tions reveal a population of bright X-ray point sources, including nine ULX, defined by 
L x > 10 39 ergs s _1 in the 0.3-8.0 keV band. The brightest source has an 0.3-8 keV lu- 
minosity L T =2xl0 40 ergs s -1 , placing it at the upper end of the X-ray binary lum inosity 



function ( jGrimm. Gilfanov. &: Sunyaevl l2003t ISwartz. Ghosh. Tennant. fc Wul 120041 ) . The 
spectrum of this source is hard, and similar to other binaries in an "ultralum inous" state 
jGladstone. Roberts, fc Donell2009l : iBerghea. Weaver. Colbert, fc Roberts! l2008h . 



We compare the NGC 922 X-ray source population with that of the Cartwheel. These 
two systems are similar in many respects (both drop-through ring galaxies) but the Cartwheel 
has a lower metallicity than NGC 922. We find that the number of ULX in NGC 922 and the 
Cartwheel scales with the star formation rate: there is no evidence for an excess of sources 
in the Cartwheel. 

Population synthesis models using StarTrack suggest that the ULX population in both 
NGC 922 and the Cartwheel is dominated by systems with strong wind accretion from su- 
pergiant donors onto direct-collapse BHs. The default models over-predict the number of 
observed ULX in both galaxies. There are several simplifying assumptions in the simula- 
tions which could account for the discrepancy. For example, the models assume that direct 
collapse events impart no natal kick to the resultant black hole. However, even small kicks 
could separate the widely spaced progenitors before the ULX phase begins, hence reducing 
the number of ULX. The simulations correctly predict the ratio of the number of ULX in 
NGC 922 and the Cartwheel. We emphasize that this relation falls naturally out of the 



Linden et al.l ( 120101 ) simulations, as the ULX population at both metallicitie s is created 
through the same evolutionary channels. We note that the iLinden et al.l ( 120101 ) simulations 
predict a substantially different ULX population below Z <0.15 Z Q . Thus it would appear 
that the metallicity of the Cartwheel is not low enough to see a difference in the ULX pop- 
ulation compared to NGC 922. We are currently working on a project to search for ULX in 
Extremely Metal Poor Galaxies (XMPG). These are the most metal poor galaxies known, 
and a logical place to find ULX if they favor metal poor systems. We will report on these 
findings in a future paper. 
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Fig. 1. — H-q image of N GC 922 with positions of X-ray sources marked. The H-a image is 
from iMeurer et al.l (120061 ). downloaded via the NASA Extragalactic Database. 
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08.0 07.0 06.0 2:25:05.0 04.0 03.0 02.0 

R.A. 

Fig. 2.— RGB Chandra image of NGC 922. The three bands are red: the 0.3-1.0keV; 
green: 1.0-2.0; and blue: 2.0-8.0. An adaptive smoothing algorithm has been applied to 
the data. The extraction regions for the point sources are defined by the white ellipses. 
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Fig. 3. — X-ray luminosity function of the Cartwheel galaxy (top, black) and NGC 922 
(bottom, red). Two fits are shown for each galaxy: the upper fit is to the raw data, an d the 
lower fit is corrected for background sources using the results of iGiacconi et al.l (120011 ) . The 
vertical lines show the lower limit to the fit for the Cartwheel (black) and NGC 922 (red). 
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